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Abstract 

We investigate complex PT and non-PT-symmetric forms of the generalized Woods- 
Saxon potential. We also look for exact solutions of the Schrodinger equation for the PT 
and/or non-PT-symmetric potentials of the kind mentioned above. Nikiforov-Uvarov 
method is used to obtain their energy eigenvalues and associated eigenfunctions. 
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1 Introduction 



x 



A large variety of potentials with the real or complex forms are encountered in various fields 
of the physics. A consistent physical theory of quantum mechanics in terms of Hermitian 
Hamiltonians is constructed on a complex Hamiltonian. In this case, its energy levels are 
real and positive as a consequence of PT-symmetry. Where P and T stand for the parity (or 
space) and time reversal operators, respectively. It is also well-known that PT-symmetry 
does not leads to completely real spectrums, because there are several potentials where part 
or all of the energy spectrums are complex. Exact solution of the Schrodiger equation for 
these potentials are generally of interest pflEHH]. Recently? Bender and his co-workers have 
studied a number of complex potentials on PT-symmetric quantum mechanics. They have 
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showed that the energy eigenvalues of the Schrodinger equation are real when PT-symmetry 
is unbroken, whereas they come in the shape of complex conjugate pairs when PT-symmetry 
is spontaneously broken |4j. In these studies, some numerical and analytical techniques have 
been used to investigate non-Hermitian Hamiltonians with real or complex spectra [SUHj- 

PT invariant operators have been analyzed for real and complex spectra by using a va- 
riety of techniques such as variational methods [7j, numerical approaches jHj, semiclassical 
estimates , fourier analysis [TU] and group theoretical approach with the Lie algebra [TT] . It 
is pointed out that PT invariant complex- valued operators may have real or complex energy 
eigenvalues. Many authors have studied on PT-symmetric and non-PT-symmetric non- 
Hermitian potential cases such as flat and step potentials with the framework of SUSYQM 
j!2j . exponential type potentials [T^j, quasi exactly solvable potentials ^3], complex Henon- 
Heiles potential [J3] and deep potential to describe optic-model analysis of elastic and in- 
elastic scattering processes [TT)] . 

Recently, an alternative method which is known as the Nikiforov-Uvarov (NU) has been 
introduced in solving the Schrodinger equation. The solution of the Schrodinger equation 
for well-known potentials and Schrodinger-like (i.e. Dirac and Klein-Gordon) equations for 
a Coulomb potential have been obtained by using this method [Tj\. The increasing interest 
on this solution method shows that PT-symmetric potentials are also suitable in solving the 
Schrodinger equation for a exponential-type potential. This type potential has been used 
in the work of Berkdemir et.al, for the nuclear scattering applications. The potential 
known as Woods-Saxon has been generalized with an additional derivative term and solved by 
using the NU method to obtain eigenvalue equations. However, PT and non PT-symmetric 
versions of the generalized Woods-Saxon potential have not been introduced in the literature 
and solved analytically for the Schrodinger equation. From this point of view, this method 
avoids direct solving of Schrodinger equation for these potential forms and makes the problem 
more interesting. 

The organization of the paper is as follows. In Sec. II, we introduced a short summariza- 
tion of the NU method, since the detailed description has been already given in Ref . JT3 and 
also used in Ref.[JS]. In Sec. Ill, the method utilized to solve the Schrodinger equation for 
PT and non-PT-symmetric non-Hermitian forms of the generalized Woods-Saxon potential. 
Finally, in Sec. IV, we concluded the article with a summary of our main remarks. 
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2 Overview of the Nikiforov-Uvarov Method 



The method of solving second-order differential equations developed by Nikiforov-Uvarov 
has been used in appearing many problems of quantum mechanics. In principle, the method 
makes it possible to present the theory of special orthogonal functions ^J] by following a 
differential equation of the form 

ns) + T ^'(s) + ^ms) = o, (i) 

a(s) a z [s) 

where a(s) and a (s) are polynomials, at most second-degree, and r (s) is a first-degree 
polynomial. Following the method used in Ref.|18j. we transform the equation for ip(s) to 
an equation of hypergeometric type 

a(s)y" + r(s)y' + Xy = 0, (2) 

by inserting ip(s) = 4>(s)y(s), where <fi(s) satisfies the equation <^(s)'/0(s) = ir(s)/a(s). 
t(s) —t (s) + 2tc(s) and its derivative has to be negative. In the present case, the polynomial 
tt(s) is 



*'(*)- r (s) ^ 

IT = ± 



\ 



2 

a'(s)— t (s 



- a (s) + ka{s) (3) 



2 

where the parameter k is a constant (k = A — 7r'). y(s) is the hypergeometric type function 
whose polynomial solutions are given by Rodrigues relation 

^l^"^ < 4 > 

where B n is a normalizing constant and the weight function p(s) must be satisfied the case 

^Hs)p(s)] = r(s)p(s). (5) 

On the other hand, the energy eigenvalues of Eq.(j3J) are determined by 

/ n(n — 1) n 

\ = \ n = —nr — -a , {n = 0, 1, 2, ...). (6) 

Hence, the solution of second-order differential equations can be solved by means of this 
method and found their energy eigenvalues and associated wavefunctions analytically. 
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3 Calculations for the Generalized Woods-Saxon Po- 
tential 



The interactions between nuclei are commonly described by using a potential that consist of 
the Coulomb and the nuclear potentials. The nuclear potential is usually taken in the Woods- 
Saxon potential form, which is one of the important potentials in nuclear physics. Moreover, 
the Woods-Saxon potential can be also used to describe the interaction of a nucleon with a 
heavy nucleus. In this paper, we selected a most general form of the Woods-Saxon potential 
given by [THj ; 



where z = exp ( x ~ a R ° ) ; Vo is the potential depth, R is the width of the potential or the 
nuclear radius and the parameter a is the thickness of a surface layer in which the potential 
falls off from V = outside to V = — V inside the nucleus. The second term of the right-side 
of Eq.(0) denotes the derivative of Woods-Saxon potential responsible for the generalization 
while the first term of the right-side of the same equation represents the Woods-Saxon 
potential. In order to calculate the energy eigenvalues and the corresponding eigenfunctions, 
the potential function given by Eq.(J7J) is substituted into the one-dimensional form of the 
Schrodinger equation: 



Here, some assignments are made in the Schrodinger equation such as \ja = 2a, q = 
exp(—2aRo) and C = 2aVoq. 

We used the NU method to obtain the exact solutions of the SE with the PT-/non-PT- 
symmetric potentials only for the s-states. The energy eigenvalues and eigenfunctions are 
found in the real or complex forms and in terms of Jacobi polynomials, respectively. 

3.1 PT-symmetric and non-Hermit ian Woods-Saxon case 

We are going to consider different forms of the generalized Woods-Saxon potential, namely at 
least one of the parameters is imaginary. For a special case, we take the potential parameters 
in Eq.(j7J) as Vq — > Vq and a —>■ iai, where aj is a real parameter of the imaginary part. Such 




(7) 
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a potential is called as PT-symmetric and also non-Hermitian, since the property V(— x)* = 
V(x) is exist. Hence, the concept of PT-symmetry can be also used in one-dimensional 
quantum mechanical problems as many problems exhibiting PT-symmetry • in this 

case, the new shape of the potential in one- dimensional space becomes 

V Ce 2iaix 



V(x) = -■ 

1 + qe 2ia i x (1 + qe ^i^f 

or it can be written as a complex function 

1+5 cos 2a jx — iq sin 2aix 



(9) 



V(x) = -V 



1 + q 2 + 2q cos 2aiX 



q ^2q + (1 + q 2 ) cos2«/X — i(q 2 — 1) sin2a/x\ 



(10) 



(1 + q 2 + 2q cos 2ajx) 

The type of this potential is known as a complex periodic potential having PT-symmetric 
and its form is given by V(x) = zsin 2n+1 (x), (n = 0,1,2,...) [2S]. A detailed discussion 
exhibits for this potential that it has real band spectra from Ref.[25]. In our case, we will 
consider the form given in Eq.(J5]l. following a procedure similar to the previous section. 
Now, in order to apply the NU- method, we rewrite Eq. (jHJ) by using a new variable of the 



form s 



2iajx 



d 2 ip(s) ldip(s) m 



ds 2 



+ 



s ds 2h 2 a 2 s 2 



E + 



Va 



Cs 



l-qs) (1-qsY 



By introducing the following dimensionless parameters 



-^>0 (P<0), p 



2h z aj y 11 ' 2h z a 2 v 11 ' 2h z aj 

which leads to a hypergeometric type equation defined in Eq.flTJ): 



if)(s) = 0. 



mC 

(7>0) 



'ID 



(12) 



d 2 ip(s) ^ 1 — qs dip(s) ^ 



1 



e(l - qs) 2 - p(l - qs) + 7s ip(s) = 0. (13) 



ds 2 s(l — qs) ds s 2 (l — qs) 2 
After the comparison of Eq. (jl3j) with Eq.Q, we obtain the corresponding polynomials as 



r(s) = l — qs, a(s) = s(l - qs), a (s) = eq 2 s 2 - {2eq - f3q - 7)5 - (3 + e. (14) 
Substituting these polynomials into Eq.Q, we obtain ir(s) function as 

tt(s) = ± - \l{q 2 - Aeq 2 - Akq) s 2 + 4 {-(3q - 7 + 2eq + k) s + 4 (J3 - e) (15) 
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taking cr'(s) = 1 — 2qs. The discriminant of the upper expression under the square root has 
to be zero. Hence, the expression becomes the square of a polynomial of first degree; 

{-(3q - 7 + 2eq + kf - ((3 - e) (q 2 - Aeq 2 - Akq) = 0. (16) 

When the required arrangements are done with respect to the constant k, its double roots 



are derived as k^ 2 = (7 — (3q) ± ((3 — e) (1 — -^). 

Thus substituting, these values for each k into Eq. ffTo]) following possible solution is 
obtained for tt(s) 



7r(s) 



(2^7 + 1 - % ) qs _ 



2 2 



for k=(<y-j3q)-qy/(p-e){l-%) 



(17) 



[ for A;=( 7 -/?g) + g v /(/3-e)(l-f) 

After appropriate choice of the polynomial 7r(s) and k, we can write the function r(s) which 
has a negative derivative as follows 



t(s) = l + 2Jp-e-qs(2 + 2J/3-e+Jl- 4 7 /g 



and then its negative derivatives become 



r'(s) = -q[2 + 2J(3-e + Jl-4 1 /q 



(19) 

A particularly interesting result of Eq. (jl9|) is that the polynomial t(s) is a generalization of 
the NU method to the complex quantum mechanics. Therefore, from Eq.fjUJ) and Eq. (jl9j) . 
we write 

\ = \ n = nq(2 + 2^(3-e + ^l- 4 7 / g ) + n(n - l)q, (20) 

and also obtain 

A = (7 - fiq) ~ qy/(J3-e)(l-4ry/q) - \ - q (yfiT^e + i^/l-47/g) • (21) 

With the comparison of A values in Eq. (j20j) and Eq. (|21|) , we found the energy eigenvalues as 
follows 



E 



nq 



2h 2 



rn 



( 



-i [l + 2n+ Jl-2mC/qh z ) + 



mV /2h 2 ai 



l + 2n + Jl- 2mC/qh 2 
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and 

g „, = f - j, (l + 2n + y/i - 4^) 2 - ( 1+2 „ + ^ ) 2 ■ (23) 

It is clear that the energy eigenvalues have a real part and then we can say that the real part 
of energy eigenvalues in Eq. (}23|) determines the energy spectrum in frame of PT-symmetric 
quantum mechanics. This situation can be also seen from the condition e nq > 0, since the 
energy eigenvalues of the generalized Woods-Saxon potential are negative. In this sense, the 
number of discrete levels is finite and determined by the inequality f|) + (j) < |, where 
b = 1 + 2n + J I — 47/g, (3 = ™Y°2 and 7 = ™f 2 . Hence, we can write a condition on 
the discrete levels for the negative energy spectrum of the complex parameter generalized 



Woods-Saxon potential if n < ymV /2h 2 a 2 I - J 1/4 - j/q - 1/2. 

Let us now find the corresponding eigenfunctions. The polynomial solutions of the hy- 
pergeometric function y(s) depend on the determination of weight function p(s) which is 
satisfies the differential equation [a(s)p(s)] f = r(s)p(s). Thus, p(s) is calculated as 

-1 Q 2JT^ 



P (s) = (1 - qs y-' swv-z, (24) 



where v = I + Jl — ^. Substituting into the Rodrigues relation given in Eq.(jU), the 



eigenfunctions are obtained in the following form 

1 r\ n 

y nq (s) = A n (i- q sr^- 1 Ks~^±- 



n+u-l n+2*Jl3- 



qs) s 



(25) 



where A n = l/n\. Choosing q — 1, the polynomial solutions of y n (s) are expressed in terms 
of Jacobi Polynomials [constant] Pn^~^' V ^(1 — 2s) which is one of the orthogonal 
polynomials. By substituting tt(s) and <r(s) into the expression 0(s)'/0(s) = 7r(s)/a(s), the 
other part of the wave function is found as 



<f>(a) = (1 - s) (26) 

where p = vjl. Combining the Jacobi polynomials with <j)(s), the s— wave complex functions 
are constructed as 

Ms) = B nS V^(l - sYP { n^ u ~ 1] (l - 2s), (27) 
where B n is the normalizing constant. 
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3.2 Non-PT symmetric and non-Hermitian Woods-Saxon case 



In order to be more specific, we are going to take the potential parameters given in Eq.fJJJ) 
as Vo — > Vor + iVoi and a — > a, where Vor and Voi are the real parameters of the complex 
part. Such a potential is called as non-PT-symmetric but non-Hermitian. In this case, the 
condition of PT-symmetry is not occurred due to the fact that V(— x)* ^ V(x). Hence, its 
shape becomes 



V{x) 



Vor Cz 



:28) 



i + z (i + zy i + z_ 

and these types potentials are called as the pseudo-Hermitian. To simplify the form of the 
above equation, we take again an independent variable by choosing the form s = —e 2ax , and 
then we obtain the generalized equation of hypergeometric type, 

1 



d 2 ip(s) 1 — qs dif)(s) 
ds 2 s(l — qs) ds s 2 (l — qs) 2 

for which 



-£(1 - qs) 2 + ((3 + iS)(l - qs) - 7 s ^(s) = 0,(29) 



r (s) = 1 — qs, cr(s) = s(l — qs), a (s) = —eq s + (2eq — (3q — idq — 7)5 + (3 + %8 — e, 

mC 



mE >o, p mVoR 



2h 2 a 2 " H 2h 2 a 2 ' 7 2h 2 a 2 ' " 2h 2 a 2 

Following the solution procedures of the NU— method, we can derive the possible solutions 
for 7r(s) as below 



<S = ^ (P, 7, 5 > 0). 



71 IS) = ± - 

v ; 2 2 



2yJe-(/3 + i6) T y 1 + y j qs ± 2yje- (/3 + i8) 



(31) 



for k = (f3q — 7 + idq) ± qJ(e — (3 — i5)(l + —). After performing an appropriate choice for 



k and 7r(s), we can write r(s) and r' as 



r(s) - 1 + 2Je - (p + iS) - qs 2 + 2Je -{(3 + iS) + 4 7 /g 



t'(s) = -q(2 + 2Je-((3 + iS) +Jl+ 4 7 /g 



(32) 



In the present case, the eigenvalue equation given in Eq.© is established as follows 



A = X n = nq f 2(1 + y/e-(P + iS) ) + J 1 + y j + n(n - l)q, (33) 
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and then if it is compared with the another form of A 
^ = (A? - 7 + iSq) - q 
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4 7 



-flV«-"-* + V 1 + 7 )• 



(34) 



the energy eigenvalues are reduced to the following form 

2 4(^) 2 (1&-1&) 



ft. 2 



1 



nq 



2ma 2 I 16 



8mCa 2 . 
1 + — o — + (l + 2n) 



1 + ^f 1 + (1 + 2ra) 



ft 2 



— / 



4(^) 2 Wo/ 



2ma? 



+ 



mVoia 2 



(35) 



1 + + (1 + 2n) 

In order to mention from an acceptable energy spectra in Eq.([35|h we have to provide the 
visible condition given by if and only if n > 8 (5 2 — (3 2 ) 1 ^ — k/2 for both Vqr > Vqj and 
V 0R < V QI , where k = yfl + ^ + 1. 

Now, our procedure is that the corresponding unnormalized wave functions should be 
determined in terms of Jacobi polynomials. By considering the Eq.(|5]) and using Eq. (J32j) . 
we obtain 

p(s) = (l-qs) R - 1 s 2 ^/ £ -P- iS . (36) 

After that, we derive the polynomial 0(s) from the equality 0(s)'/0(s) = 7r(s)/a(s) as 
follows 

(f)(3) = s ^ £ - p - iS (l - qs) K/2 . (37) 
Hence, we found the relevant wave functions in terms of Jacobi polynomials 

Ms) = D nS V^\i - s) «/2p,?V^^ «-d (1 _ 2s); (38) 

where D n is the another normalizing constant and the parameter q is chosen as 1. 

4 Conclusions 

The exact solutions of the Schrodinger equation for the generalized Woods-Saxon potential 
with zero angular momentum are obtained to found the energy eigenvalues and eigenfunc- 
tions. Nikiforov-Uvarov method is used to solve the Schrodinger equation. The correspond- 
ing wave functions are expressed in terms of Jacobi polynomials. It is pointed out that 



9 



the NU method can be generalized by the complex polynomials n(s) and r(s) with the 
framework of PT-symmetric quantum mechanics. When the potential parameter a is purely 
complex, it is seen that the energy eigenvalues are in the real form. In this point, the shape 
of potential takes the complex PT-symmetric form. It is not difficult to see that the en- 
ergy eigenvalues determines the discrete levels or energy spectrum under the condition of 
n < \J mVo / 2h 2 a] — ^/l/4 — 7/g — 1/2. Therefore, if the parameter of potential Vo is chosen 
as complex, it is clear that the energy eigenvalues represent the complex energy spectrum. 
We realize that the energy eigenvalues consist of essentially a real component. Consequently, 
we find that the energy eigenvalues of a single particle within the complex PT and non-PT- 
symmetric potentials are given in Eq.([22p and Eq. (j35|) . respectively. 
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